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Abstract. We use the three-dimensional Heisenberg model with site randomness as an 
effective model of the compound Sr(Fei_ a; Mn a; )02. The model consists of two types of ions 
that correspond to Fe and Mn ions. The nearest-neighbor interactions in the afe-plane are 
antiferromagnetic. The nearest-neighbor interactions along the c-axis between Fe ions are 
assumed to be antiferromagnetic, whereas other interactions are assumed to be ferromagnetic. 
From Monte Carlo simulations, we confirm the existence of the double-q ordered phase 
characterized by two wave numbers, (7T7r7r) and (7T7tO). We also identify the spin ordering 
pattern in the double-q ordered phase. 



1. Introduction 

Infinite-layer iron oxide SrFe02 exhibits a Neel transition to the G-type antiferromagnetic 
ordered phase at T/v = 473K[U [2]. The transition temperature decreases as Mn 2+ ions 
are substituted for Fe 2+ ions[3]. In the compound Sr(Fei_ a; Mn a; )02, Fe and Mn ions are 
magnetic and placed randomly on a tetragonal lattice. In neutron scattering measurements on 
Sr(Feo.7Mno.3)02, magnetic peaks have been observed to develop at two distinct wave vectors 
q = (mm) and (7T7tO). These orders coexist at low temperatures |3j. However, the spin ordering 
pattern in the low-temperature phase is not known. 

The site-random model has been used as an effective model of the random magnets. The 
site-random model consists of two types of ions labeled A and B. The ions are placed randomly 
on the lattice. The interactions of each bond are set depending on the combination of ions. In 
previous works that employed the site-random model [U EJ [TJ [8], all interactions depended on 
the arrangement of A and B ions: the interactions between A ions were ferromagnetic, whereas 
other interactions were antiferromagnetic. In the case of the Ising spin model on a simple cubic 
lattice[4J, it was found through Monte Carlo simulations that the phase diagram lacks both a 
mixed ordered phase and spin-glass phase. In the Heisenberg spin model[5l[6j[71[8], as m ^ ne 
spin case, there is also no spin-glass phase in the phase diagram. However, the ferromagnetic 
and antiferromagnetic mixed ordered phase exists at low temperatures. In this mixed ordered 
phase, the peaks of the structure factor develop at wave vectors q = (000) and (-ett-k). Moreover, 
the ferromagnetic and antiferromagnetic ordering vectors are mutually perpendicular. For this 
model, the occurrence of a two-step phase transition has been reported: one is the transition 
from the paramagnetic phase to the ferromagnetic (antiferromagnetic) ordered phase, the other 
is the transition from the ferromagnetic (antiferromagnetic) ordered phase to the mixed ordered 
phase. The universality class at each phase transition has been suggested to belong to the 
three-dimensional Heisenberg universality class [8j. 



In the previous works that used the site-random Heisenberg model [5j El [8], q = (000) and 
(7T7T7r) orderings have been found. However, since the magnetic peak positions are q = (tttttt) and 
(7T7rO) in Sr(Feo.7Mno.3)02, the spin ordering pattern is different from that observed previously. 
Thus, our current aim is to clarify the spin ordering pattern in the low-temperature phase of 
Sr(Fei_ a; Mn ;c )0 2 . 

2. Model 

To introduce suitable randomness into Sr(Fei_ :r Mn 2 ;)02, we investigate a classical Heisenberg 
model with site-random interlayer couplings on a simple cubic lattice: 

H = -J ^2 Si-Sj- ^2 JijSi-Sj, (\Jij\ = \J\)- (1) 

-plane (id) 

The first term of the Hamiltonian denotes the uniform antiferromagnetic nearest-neighbor 
interaction (J < 0) in the a6-plane. The sign of the nearest-neighbor coupling Jij in the second 
term depends on the arrangement of two ions along the c-axis; the interactions between Fe 
ions are assumed to be antiferromagnetic (Jjj < 0), whereas other interactions are assumed to 
be ferromagnetic (Jy > 0). For simplicity, we further assume that the absolute values of all 
interactions are the same and that the spin lengths of Fe and Mn ions are the same. Hereinafter, 
we set | J | = 1 and denote the concentration of Mn ions by pMn- From the above assumptions 
regarding the interactions, at pyin = where all magnetic ions are Fe ions, the system exhibits 
a Neel transition to the G-type antiferromagnetic ordered phase characterized by q = (tttttt), 
because all interactions along the c-axis are antiferromagnetic. In contrast, at pyin = 1 where 
all magnetic ions are Mn ions, the C-type antiferromagnetic ordered phase characterized by 
q = (7T7rO) becomes stable at low temperatures, because all interactions along the c-axis are 
ferromagnetic. For < pyin < 1> the frustration is caused by the random ferromagnetic and 
antiferromagnetic interactions along the c-axis direction. 

3. Simulation Results 

We use Monte Carlo simulations based on the standard heat-bath method to the model on 
&nN = LxLxL simple cubic lattice with a periodic boundary condition. Before starting 
the simulations, we set pyin so that the number of Mn ions is an integer, and we set random 
configurations of Fe and Mn ions according to pMn- In this paper, pMn is set to 0.1875. For this 
value of pyin, there are more antiferromagnetic couplings than ferromagnetic couplings along the 
c-axis. We prepare 64 samples, and each run contains 10 6 -10 Monte Carlo steps per spin at 
each temperature. 

The temperature dependence of the structure factor S(q) at q = (7r7T7r) and (7T7tO) is shown 
in Figure [TJ The structure factor S(q) is defined by 

where (• • • ) indicates the thermal average. 5(7r7T7r) and S(tttt0) increase at T = 1.09 and 0.46, 
respectively. From Figure [TJ it is clear that there are three types of phases in this model. In the 
low-temperature phase (phase I), the peaks of the structure factor develop at two distinct wave 
vectors q = (tttttt) and (7T7tO). Accordingly, this phase is called a double-q ordered phase. In the 
intermediate-temperature phase (phase II), since the structure factor has a peak at q = (tttttt), 
this phase is the G-type antiferromagnetic ordered phase. The high-temperature phase (phase 
III) is the paramagnetic phase. The transition temperature from the G-type antiferromagnetic 
ordered phase to the paramagnetic phase is lower than the Neel temperature of the pure system 
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Figure 1. Temperature dependence of the structure factor at q = (iririr) and (mrO) when the 
lattice size is L = 12 and 24. Phases denoted by I, II, and III are the double-g ordered phase, 
the G-type antiferromagnetic ordered phase, and the paramagnetic phase, respectively. 

T/v — 1.44 [9j . From the results of the structure factor, we conclude that the double-g ordered 
phase characterized by the wave vectors q = {tttttt) and (7T7rO) exists in the model. The positions 
of the magnetic peaks for the model are the same as those observed for Sr(Fei_ x Mn x )02- 

@ 

To identify the spin ordering pattern in the double-g ordered phase, we calculate the 
correlation function along the c-axis between two spins. The correlation function along the 
c-axis is given by 

G c (r 2 ) = ^( S (r i )-s(r i + r z e 2 )), (3) 

i 

where r z is the distance between two spins, and e z is unit vector of the c-axis direction. Figure [2] 
shows dependence of the correlation function on distance r z for L = 24 at T = 0.1. Since the spin 
correlation is ferromagnetic for an even number of r z , all spins in odd (even)-numbered layers 
are co-linear; The Neel-like configuration of each layer is stacked in parallel to other layers. 
Furthermore, the spin correlation is weakly coupled antiferromagnetic for an odd- number of r z . 
This suggests that the angle between nearest-neighbor spins along the c-axis is close to 2ir/3 at 
T = 0.1. 

To investigate the temperature dependence of spin ordering pattern, we calculate the angle 
6 between the staggered magnetization of an odd-numbered layer, M , and that of an even- 
numbered layer, M e . M Q and M e are given by 

M = | £ e^s(r), (4) 

rSodd-numbered 
layer 

M e = A Y, e iq - r *(r), (5) 

rGeven-numbered 
layer 

where q = (-7T7rO). We define the angle 9 between M Q and M e as follows: 
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Figure 2. Correlation function along the c-axis for L = 24 at T = 0.1. 



CD 



1 

0.95 
0.9 

0.85 
0.8 

0.75 
0.7 



L=12 — - 
L=16 — b- 
L=20 — •- 
L=24 — »- 




fita • 



II 



III 



■ 
• 



0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 3. Temperature dependence of the angle 8 between the staggered magnetization of an 
odd-numbered layer and that of an even-numbered layer. 



Figure [3] shows the temperature dependence of 8. At the higher-temperature transition point, 
the value o£8/ir approaches 1 and the (7T7T7r) order develops. At the lower-temperature transition 
point, the value of 8/n deviates from 1. In phase I, the value of 8 depends on the temperature 
and monotonically decreases as the temperature decreases. Thus, the spin configuration 
characterized by 8 changes depending on the temperature in phase I. At a temperature of 
zero, the value of 8/tt appears to converge to 0.73. From the results, we find that the lower- 
temperature phase transition is caused by 8/tt deviating from 1. 

Let us consider the translational symmetry and the spin rotation symmetry in each phase. 
The intermediate-temperature phase is translationally symmetric with the 0(3) spin rotation 
symmetry broken down to the U(l). In the double-g ordered phase, both the translational 
symmetry and the U(l) spin rotation symmetry are broken. Thus, we conclude that the 
transition to the double-g ordered phase is characterized by the breaking of the U(l) spin 
rotation symmetry. 



4. Discussion 

We discuss dependence of physical properties on the concentration of Mn ions pMn- In this 
paper, we have fixed pMn = 0.1875. In this case, there are more antiferromagnetic couplings 
than ferromagnetic couplings along the c-axis, and the intermediate G-type antiferromagnetic 
ordered phase exists in a broad temperature range (0.46 < T < 1.09). At pyin — 0.3, the number 
of antiferromagnetic couplings and ferromagnetic couplings along the c-axis become nearly the 
same. Then, (ttit-k) order and (ttttO) order develop at nearly the same temperature, and the 
intermediate-temperature phase exists in a very narrow temperature range. Furthermore, for 
PMn > 0.3, the intermediate C-type antiferromagnetic ordered phase exists, and the transition 
to the double-g ordered phase occurs at low temperature. In future work, we will attempt to 
clarify the details of the phase diagram of the temperature versus the concentration of Mn ions. 

5. Conclusion 

We have used the three-dimensional Heisenberg model with site-random interlayer couplings as 
an effective model of the compound Sr(Fei_ x Mn x )02. The model consists of two types of ions 
that correspond to Fe and Mn ions. The interactions have been assumed to be antiferromagnetic 
in the a6-plane. The interactions along the c-axis between Fe ions have been assumed to be 
antiferromagnetic, whereas other interactions have been assumed to be ferromagnetic. From 
Monte Carlo simulations, we have confirmed the existence of the double-q ordered phase 
characterized by the wave vectors q = (tttttt) and (mrO). The positions of the magnetic peaks 
for the model are the same as those observed for Sr(Fei_ x Mn x )02- We have also identified the 
spin ordering pattern in the double-g ordered phase. In the double-q ordered phase, all spins 
in odd (even)-numbered layers are co-linear. Furthermore, the staggered magnetization of an 
odd-numbered layer and that of an even-numbered layer are not parallel, and the angle between 
them changes depending on the temperature. We conclude that the transition to the double-g 
ordered phase is characterized by the breaking of the U(l) spin rotation symmetry. 
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